Shock waves in two-dimensional granular flow: effects of rough walls and 

polydispersity 
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We have studied the two-dimensional flow of balls in a small angle funnel, when either the side 
walls are rough or the balls are polydisperse. As in earlier work on monodisperse flows in smooth 
funnels, we observe the formation of kinematic shock waves/density waves. We find that for rough 
walls the flows are more disordered than for smooth walls and that shock waves generally propagate 
more slowly. For rough wall funnel flow, we show that the shock velocity and frequency obey simple 
scaling laws. These scaling laws are consistent with those found for smooth wall flow, but here they 
are cleaner since there are fewer packing-site effects and we study a wider range of parameters. For 
pipe flow (parallel side walls), rough walls support many shock waves, while smooth walls exhibit 
fewer or no shock waves. For funnel flows of balls with varying sizes, we find that flows with weak 
polydispersity behave qualitatively similar to monodisperse flows. For strong polydispersity, scaling 
breaks down and the shock waves consist of extended areas where the funnel is blocked completely. 

PACS numbers: 45.70.Mg, 45.70.Vn 



I. INTRODUCTION 

Density waves can occur in various granular flow sys- 
tems such as funnels and hoppers jl) , pipes |2) , and hour 
glasses ra. In previous experiments, we studied the for- 
mation of kinematic shock waves which propagate against 
the main flow in a two-dimensional system of balls rolling 
in a small angle funnel j|, ||, |(| . A sketch of the setup used 
is shown in Fig. |l|(a)-(b). In this earlier work, the walls of 
the funnel were smooth, and the balls were of equal size 
(monodisperse flow) . In the present work, we study what 
happens when we break the peculiarities of this smooth 
wall/monodisperse system by either (i) making the walls 
of the funnel rough (Fig. |l](c)), or (ii) taking "polydis- 
perse" mixtures of balls of different sizes (Fig. [l](d)). 

The crucial experimental parameters characterizing 
the flow geometry are the funnel opening angle (5 and the 
funnel outlet width D. For smooth wall/monodisperse 
flows the most important features of the shock waves 
were found to be the following g §, g|: (i) For (3 > 0, 
the rolling grains tend to locally form triangular lattices 
which lead to the creation of shock waves predominantly 
at particular sites in the funnel where close packing oc- 
curs, (ii) The velocities U of the shocks are, in good 
approximation, a function of the rescaled width w{x)/D 
only. Here x is the coordinate along the funnel (Fig. 0(b)) 
and w(x) denotes the funnel width at position x. 

As we will show below, by making the walls rough or 
the flow polydisperse, the triangular packing can (par- 
tially) be suppressed. We have studied the shock statis- 
tics and the behavior of individual balls as in previous 
work ||, ||] , and will present here the results for the cre- 
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FIG. 1: (a) Sketch of the experimental setup, (b) Schematic 
top view of the setup, showing the important parameters of 
the funnel geometry, (c) Close-up of the balls rolling between 
rough walls. The linked balls near the edges are glued to the 
walls (the scale is in cm), (d) Snap shots of two different 
polydisperse mixtures, referred to as mix I and IV (see text). 



ation and propagation of shocks in these systems in the 
pipe flow (parallel walls) and in the small angle, inter- 
mittent flow regime ||. We have also investigated, us- 
ing ball tracking, the formation of shear bands for rough 
walls and the effects of completely stationary shock pack- 
ings in polydisperse flows. 

Recently, a number of related experiments on two- 
dimensional flow have been performed. Tsai et al.Q| 
studied a system close to our rough wall pipe flow exper- 
iment, and included the effect of partially blocking the 
outlet. Reydellet et al. || studied a falling vertical col- 
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umn of balls where ball-ball interactions are dominated 
by collisions, and the rolling of balls does play a minor 
role. Finally, Le Pennec et al. |)| studied two-dimensional 
rolling flows of small glass balls in flow geometries with 
very large funnel angles. 

The experimental setup and methods of analysis have 
been described in detail elsewhere [0, |6| and we summa- 
rize the essential aspects here. A sketch of the setup is 
shown in Fig. [j](a) and the important geometric param- 
eters are shown in Fig. |(b). The balls roll on a coated 
Lexan plane (with inclination 9 = 4.1°) in a single layer, 
in a funnel formed by aluminum walls and covered by a 
transparent lid. The granular material is comprised of 
50000 brass balls of 3.18 mm diameter in the monodis- 
perse case, while details of the polydisperse mixtures are 
discussed in section [II. The rough walls were made by 
glueing linked rows of balls of nearly the same diameter as 
the rolling balls to the original smooth walls (Fig. 0(c)). 
The walls are straight for 200 cm and curve smoothly at 
the top to form a reservoir. They can be moved to vary 
the outlet width D (0-35 mm) and the funnel half angle 
(3 (0° — 3°). A light box is placed below the funnel to 
illuminate the balls from below, and a video camera is 
placed above the system. Snapshots of a small part of 
the funnel show the effect of rough walls (Fig. 0(c)) and 
polydispersity (Fig. 0(d)) on the packing of the balls. 

The paper is structured as follows: In Section || we 
present data based on monodisperse flows between rough 
walls. Shock wave statistics for both intermittent flows 
((3 > 0°) and for pipe flows (/3 = 0°) are discussed, and we 
also study data based on ball tracking. In Section III we 



investigate the effect of various degrees of polydispersity 
in flow between smooth walls. Finally, in Appendix , we 
discuss some similarities with traffic flow. 



II. ROUGH WALLS 

The basic phenomenology of the formation of shock 
waves is illustrated in Fig. 0, which shows ten subsequent 
snapshots of the balls in a small section near the outlet 
of the rough wall funnel. In dense regions, which may oc- 
cur due to a combination of geometric effects (the finite 
funnel angle (3) and the inelastic nature of the collisions, 
kinetic energy is dissipated rapidly and so the balls here 
have a lower velocity than the balls in dilute regions. This 
leads to kinematic shock waves, in which the balls become 
almost stationary and tend to pack in a lattice which ex- 
tends from wall to wall. In such a region, both energy 
and momentum are efficiently transferred to the walls. 
Shocks grow in the upstream direction due to incoming, 
high velocity balls, and dissolve in the downstream di- 
rection due to the action of the (effective) gravity which 
eventually accelerates the balls again. To analyze these 
shocks, we recorded ~ 2 — 9 second movies of 37 cm sec- 
tions of the flow (usually at 220 frames/s), and applied 
ball tracking software [||, gaining detailed knowledge of 
the position, velocity and acceleration of the individual 
balls (see Section II C ) . 

In many cases, however, one may be more interested 
in the overall features of the shocks and not so much in 
the individual balls. For such analysis, sequences of im- 
ages taken at 60 frames/s and covering 100 cm sections 
of the funnel were averaged in the transverse direction to 
obtain a one-dimensional relative density profile. An ex- 
ample of a space-time plot of this relative density p(x, t) is 
shown in Fig. 0for two different sets of parameters. From 
these the position of a shock wave and its creation site, 
average local velocity U(x) and frequency v{x) can be 
determined. As with smooth walls, shock waves between 
rough walls are easily distinguished by eye || . They are 
created at various positions in the funnel and there are 
some noticeable interactions between shocks. As in ear- 
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FIG. 2: Film sequence showing 18 cm of the funnel showing 
a propagating shock wave at /3 = 0.2°, D — 20 mm (each ver- 
tical 'stripe' is a separate picture). Such movies are recorded 
at 220 fra mes/ s and are used for the ball tracking discussed 
in section [I C For clarity, only every eight frame is shown. 
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FIG. 3: Space-time plots of the density p~(x,t) at /? = 0.4° 
for, (a) D = 15 mm and (b) D = 25 mm. 
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lier work on smooth walls, shock waves are created more 
often at large funnel angles and for small values of D. 
The shock velocity U is generally observed to grow with 
increasing /?, x and decreasing D. 



A. (3 > 

We will study now in detail the statistics of shock 
creation, velocity and frequency, based on density data 
taken for D ranging from 15 to 30 mm, and for (3 rang- 
ing from 0.1° up to 1.2°. The case of pipe flow ((3 — 0), 
where rough and smooth walls s how very different be- 
havior, will be studied in section [IB, and some results 



following from ball tracking are discussed in section II C . 



(this aligns possible packing sites). Clearly, packing ef- 
fects are present for D = 20 mm, but are washed out for 
D — 25 mm. This is consisted with the D = 15 mm and 
D = 30 mm data sets, although in these cases it is dif- 
ficult to obtain a precise estimate for the shock creation 
rates. A comparison of the data shown in Fig. ^ to data 
on the shock creation rates in smooth funnels, as shown 
in the top of Fig. [l^, confirms that rough walls suppress 
the effects of packing. 

We therefore have shown that packing effects are sup- 
pressed in sufficiently wide rough wall flow, and that in 
this case shock waves are created everywhere in the fun- 
nel with equal probability. 



2. Shock velocity 



1. Shock creation 

The monodispersity of the balls permits them to form 
close packed triangular lattices at certain packing sites 
x = Xi i n the funnel. For smooth walls, the positions of 
these packing sites are given by 



Xi = 



2r + V3r(i - 1) - D 
2tan/3 



(1) 



where i is an integer and r is the ball radius ||. In 
a funnel with smooth walls, packing effects are rather 
pronounced as shown in Fig. 18 of [pi. 

For rough walls, it is slightly more difficult to give an 
estimate of where one expects packing sites. It is, for 
example, not obvious whether the value of D (which is 
the minimum distance between the two rough walls) is 
the relevant parameter. We have found that packing ef- 
fects persist from small D, but are washed out for larger 
funnel outlets. This is illustrated in Fig. |J which shows 
the manually counted shock creation rates for two dif- 
ferent values of D as a function of the local width w 
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FIG. 4: Manually counted shock creation rates (SCR) 
for rough walls shown as a function of the local funnel width 
w(x) for D = 20 mm and D — 25 mm (both data set cover 
x = - 100 cm). 



From the density fields, we can obtain the local shock 
velocity U(x) using the shock detection algorithm RGE 
(.Relative density contrast Gradient fiige detection) [ pT[ . 
This is a refinement of the GE-method described in detail 
in @], and is used to measure U(x) and v(x). The shock 
velocities obtained with this method for forty different 
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FIG. 5: Average shock velocity U(x) for 40 different sets of 
parameters /3 and D as function of (a) the local width w(x) 
and (b) the rescaled, non dimensionalized width w(x)/D. 
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sets of parameters are shown in Fig. ||(a) as a function 
of the local funnel width w(x). This already indicates a 
clustering of the data in groups given by the funnel outlet 
D. When the data is replotted against the dimensionless 
parameter w(x)/D as in Fig. ||(b), there is a fairly good 
data collapse with U(x) » 79[w(x) / D — 0.73] cm/s. A 
similar scaling (U(x) w 6A\w(x)/D — 0.14] cm/s) was 
observed for smooth walls |], q|. Experiments with one 
rough and one smooth wall, which we will refer to as 
"semi rough" (not shown) gave a more disordered flow 
with a decent data collapse of U(x) as 83[w(x)/D — 
0.57] cm/s, which falls between rough wall and smooth 
wall data for low w{x)/D (see Fig. [7(a)). Therefore we 
conclude that the roughness of the walls leads, in general, 
to a slowing down of the shock waves. 

In the data collapse Fig. ||(b) there seems to be a trend 
for wide funnels to have slightly larger velocities. Some 
of the D = 30 mm data sets had problems with static 
buildup at the outlet during parts of the experiment 
which is likely to have affected shock statistics (increas- 
ing U{x) and v{x)) near the outlet. It is also possible 
that for larger outlet widths the data for rough walls 
becomes more comparable to those of smooth wall sys- 
tems. Nevertheless, we find that the shock velocity is 
in good approximation a function of the local rescaled 
width w(x)/D only, both for smooth and rough walls. 



3. Shock frequency 

From the density plots we can also obtain the local 
shock frequencies v{x) using the RGE algorithm Q. In 
Fig. ||(a) we show the shock frequencies for our rough 
wall data as a function of w(x). While some packing site 
periodicity is still visible for D = 15 mm, for larger values 
of D the curves are fairly smooth and we will discuss the 
data-collapse that occurs there. 

From the data shown in Fig. ro(a) it appears that v goes 
to zero when w(x) — > D, i.e., near the outlet. Rcplotting 
the data then as function ofw(x)/D — l (not shown) con- 
firmed this, and indicated that to achieve data collapse, 
one needs to plot the product v(x)D as a function of 
w(x) /D—l. Finally, in contrast to the shock velocity, the 
frequency clearly grows non-linearly with w(x)/D—l. On 
a log-log plot, it appears that vD scales as (w(x)/D— l) p , 
with an exponent p around 0.7 (Fig. ^(b)). While the 
scaling range is to small to determine whether such pow- 
erlaw scaling holds asymptotically it is a useful way of 
collapsing our data, as also is shown in Fig. ||(c) Note 
that significant deviations can be seen for larger values 
of w(x)/D. The data falling below the curve in Fig. ^(c) 
mainly belong to data sets with D = 15 mm (3 > 0.5° 
where the RGE method tends to underestimate v{x) for 
large x,/3 @. 

In a previous paper Q , which discussed the shock dy- 
namics for smooth funnels, the issue of v(x) data collapse 
was not discussed [[TJ) . Reprocessing of this smooth wall 
data employing the RGE algorithm does show a decent 
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FIG. 6: (a) Average shock frequency v{x) vs. w(x). (b) 
Rescaled shock frequency v(x)D vs. (w(x)/D) ' 7 (lin/lin). 
All for 40 different sets of parameters /3 and D. 



data collapse where v(x)D = 1.9(1) (w(x)/D - l) a7 W, 
but the number of old data sets with D ^ 10 mm that 
gives reliable ^(x)-data is insufficient to judge whether 
v(x)D or v{x) gives a better data collapse. We will dis- 
cuss this issue below for wea kly p olydisperse mixtures 
in smooth funnels in Section III 2, where more data is 
available. 



Note that the shock frequency for smooth walls is gen- 
erally less than half the frequency for rough walls. How- 
ever, since we believe that shock waves are, in general, 
produced by disturbances in the flow, this should not be 
a big surprise: rough walls produce more disturbances. 

Data obtained from flows with one rough wall and one 
smooth wall ("semi rough") confirm this intuitive pic- 
ture, since we find a decent data collapse of u(x)D w 
3A(l)(w(x)/D- l) - 7 ^) cm/s (not shown). This is right 
between the fits for the rough wall and smooth wall data 
discussed above. 
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FIG. 7: Linear fits for 1 < w(x)/D < 2.5 of U(x) vs. w(x)/D 
in (a) and for v(x)D vs. w(x)/D in (b). Both graphs show 
the lines for smooth wall, smooth/rough combined walls and 
rough wall flows. Data for weakly polydisperse, smooth wall 
flows are also included. 
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FIG. 8: Velocity (a) and rescaled frequency (b) of shock waves 
in rough wall pipe flow. 



4- Rough wall intermittent flow: conclusion 

The rough walls allow us to suppress packing effects, 
so that two simple scaling laws for the shock velocity and 
frequency emerge: 



U(x) « a 1 {w{x)/D-a 2 ) 
v(x) « D- 1 a 3 (w(x)/D- 



\0.7 



(2) 
(3) 



The fitting coefficients ai,a2 and vary with the rough- 
ness of the walls. In Fig. [7] we show the best fits for data 
obtained for rough walls, semi-rough walls, smooth walls 
and smooth walls / weakly polydisperse balls (discussed 
in detail in Section |III 2| ). The fits are all made in the 
interval 1 < w(x)/D < 2.5 for consistency. 

As shown in Fig. |^(a) , changes in the boundaries have 
a weak effect on ax, but substantially affect a%. The old 
smooth wall data for the velocity seems to deviate a lit- 
tle from this trend, and we do not know the reason for 
this. For the frequency, the roughness of the walls has 
a profound influence on the prefactor a 3 . For small an- 
gle funnels, rough walls promote the occurrence of shock 
waves, but do not affect their propagation velocity sub- 
stantially. 



B. Pipe flow 

The behavior of shocks in pipe flow (f3 = 0°) can be 
expected to be rather different from those in intermittent 
flow (/3 > 0.1°) since for pipe flow there is no "geomet- 
ric" source for the formation of shock waves, and dissi- 
pative collisions are now presumably the dominant origin 
of shock wave formation. 

Indeed, for smooth walls the (few) shock waves that 
are generated never reach the reservoir and consequently 
the flow rate (and indirectly all other flow properties) are 
entirely determined by the reservoir outflow (J, 0. Sur- 
prisingly, the qualitative features of pipeflow with rough 



walls are similar to the intermittent flows. For example, 
many shocks do travel upwards and reach the reservoir. 

Similar to smooth wall pipe flow, the rough wall pipe 
flow is extremely sensitive to exact experimental condi- 
tions and consequently the statistics obtained for these 
flows is bound to be more noisy. In fact, as we will see 
below, subtle problems near the outflow area render part 
of the data unusable. The analysis of shock properties is 
further hampered by the fact that the RGE shock iden- 
tification method |1J] used above has problems detecting 
some shocks with U < 15 — 20 cm/s. We have therefore 
decided to use an alternative method, based on space- 
time correlation functions, to obtain a measure of the 
velocity, and determine their frequency based on a sim- 
ple threshold algorithm (see Appendix B). 

Based on these methods, we have determined the shock 
velocities for rough wall pipe flow. We find that the data 
for the lower part of the funnel (0 < x < 100 cm) does not 
show any systematic trend, and often here data sets show 
a complicated mixture of periods of stationary shocks and 
moving shocks. We think that this is due to experimen- 
tal problems (in particular static charge buildup near the 
funnel outlet in some measurements). In later runs, when 
we studied the upper funnel (100 < x < 200 cm), these 
problems were solved. Here we find that the shock ve- 
locity is essentially independent of the width D, apart 
from some effects near the outflow of the reservoir (see 
Fig. 1(a)). 

The shock frequencies decrease with D, but do not 
show a clear trend with x. Rescaling the frequency with 
D, we find that vD is fairly constant, for D ranging from 
15 mm to 30 mm (see Fig. ||(b)). Again, most data for 
(0 < x < 100 cm) does not show any systematic behavior, 
due to the large amount of stationary shocks. 

Using the linear fits of U(x) and v(x)D sinilar to 
Section II A, we find U{w(x) = D) « 21 cm/s and 
v(w(x) = D)D 0.8 cm/s. Considering the systematic 
differences between RGE based and space/time correla- 
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FIG. 9: Average v x as function of y in various pipe flows (a) 
and 13 > 0.1°-fiows (b). Histograms of ball positions (consid- 
ering the width of the balls) grouped according to the v^-value 
for a pipe flow (/3 = 0.0°, D — 20 mm) in (c) and intermit- 
tent flow (/3 = 0.4°, D = 20 mm) in (d). The average square 
transverse velocity < v'y >= T y is shown as function of y for 
three pipe flows in (e) and for /3 > 0° flows in (f). Note that 
in (b),(d) and (f) the rescaled coordinate y/w(x) is used on 
the horizontal axis. 



tion based measurements of U(x) there is a fairly good 
correspondence between the pipe flow shock wave data 
from the x > 100 cm part of the funnel and the shock 
wave behavior for f3 > 0° flows. This suggest that the 
mechanisms of shock wave propagation are not too dif- 
ferent after all. 

We conclude that in good approximation, rough wall 
shock flow is dominated by shocks that are created 
mainly near the outflow region, travel upward with fairly 
constant velocity and whose frequency decreases as 1/D. 



C. Ball tracking 

We have processed detailed films such as shown in 
Fig. H to obtain ball trajectories (details of this method 
are described in ||). From the full set of ball trajectories, 
we have constructed continuous one-dimensional Eule- 
rian fields of the relative density p(x,t), velocity v(x,t), 
and acceleration a(x, t)M. A comparison of the accel- 
eration fields for rough and smooth walls indicates that 



the flow for smooth walls has more disturbances and the 
shocks are less sharply defined. 

We have previously shown (Appendix C in H) that 
smooth wall flows are reasonably one-dimensional, i.e., 
the density, velocity and acceleration fields do not show 
a strong dependence on y, the coordinate across the fun- 
nel. This is no longer true for rough wall flows as shown 
in Fig. ^ The velocity < v x > as a function of the 
transverse coordinate y is shown for various pipe flows in 
Fig. ||(a) and for (3 > 0° flows in Fig. ||(b) . In both cases 
do we find that < v x > drops off near the boundaries, 
and this effect is most pronounced for pipe flow and small 
angle flow. A number of different flow behavior could 
underly these statistics, the most obvious being: (i) The 
flow has a shear component near the walls, i.e., balls near 
walls move typically slower than in the bulk, (ii) In shock 
packings the flow reaches from wall to wall, but in fast 
regions (between shocks) balls near walls are repulsed. 
In this interpretation, balls in fast regions are repelled 
from the boundaries, and this transverse momentum is 
absorbed via ball-ball interactions in the funnel center, 
leaving the regions near the walls relatively empty. 

To resolve this ambiguity we plot histograms of the y 
positions of balls with their v x in a certain interval in 
Fig. |(c) (for pipe flow) and in Fig. |(d) (for (3 = 0.4° 
flow). If reason (ii) would dominate, the histogram of the 
slowest balls (in shock regions) should be fairly flat, but 
this is not the case and therefore we conclude that there 
is an important shear component to the flow. 

The y dependence of the y component of the granular 
temperature, T y :=< Vy >, is shown in Fig. ||(e) for pipe 
flow and in Fig. §(f) for (3 > 0° flow. For the D = 15 
mm pipe flow, T y (y) is constant while the wider pipe 
flows exhibit a more quiet region in the center - perhaps 
consistent with the "transverse momentum sink" men- 
tioned in (ii) above. The (3 > 0° flows in Fig. ||(f) are 
denser and slower, showing that in dense flows ((3 = 1°, 
(3 = 2°) transverse momentum is absorbed immediately 
and the profile is flat. The "transverse heat sink" does 
not play a role here. 

III. POLYDISPERSE FLOWS 

To suppress the close packing effects that occur for 
monodisperse balls rolling in smooth funnels, we have 
also explored flows of balls of mixed sizes. We have stud- 



Mixture 


2.5 mm 


3.0 mm 


3.2 mm 


3.5 mm 4.0 mm 


I 




10000 


10000 




II 




10000 


10000 


5000 


III 


10000 


10000 


10000 




IV 




10000 


10000 


3000 



TABLE I: Quantities of balls of various sizes in the four dif- 
ferent polydisperse mixtures used here. 
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FIG. 10: Space-time diagrams of the density p(x, i) showing 
differences in shock behavior between weak and strong poly- 
dispersity at @ = 0.5°, D = 10 mm for Mix I (a) and Mix IV 
(b). 



icd shock waves for four different mixtures, which we will 
refer to as mixture I-IV in increasing order of "polydis- 
persity" (see table 1). 

A crucial and unexpected effect of polydispersity can 
be seen in the space-time plots of the density such as 
Fig. JIG. For a weakly polydisperse mixture (mix I in 
Fig. |lO|(a)) the shocks appear to display similar behavior 
as for monodisperse balls, but for increasingly polydis- 
perse mixtures, shock waves lead to the blocking of finite 
fractions of the funnel. Here all balls are stationary for a 
finite time interval. Some clear examples of this "freez- 
ing" can be see for mix IV in Fig. [l(](b), for x w 80 cm 
and t w 3.5 s. The occurrence of finite blocked fractions 
can be interpreted in terms of a competition between the 
leading and trailing edge of a shock. For strongly poly- 
disperse mixtures, the velocity of the trailing edge of the 
shock where the shocks dissolve, is lower than that of 
the leading edge where the shock grows due to incoming 
balls. As a function of its lifetime, a shock will there- 
fore spread out, and finite regions of the funnel will be 
blocked. In contrast, when one would generate an ex- 
tended blocked area in a monodisperse flows, it appears 
that the velocity of the trailing edge is larger than that of 
the leading edge. Such a shock would then shrink, until 
leading and trailing edge come very close together and 
the shock looses its spatial extent. 

We have measured, employing ball tracking methods, 
the velocities of individual balls in shocks, and some rep- 

(a). The ball 
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resentative results are shown in Fig. 
tracking method has not been specifically fine tuned for 
polydisperse flow, which leads to a slightly noisier deter- 
mination of v x . The overall picture that emerges is that 
in monodisperse flows, balls almost never stop completely 
in a shock (notice that for the example in Fig. [ll](a) the 
minimum velocity of the balls in monodisperse shocks in- 
deed stays finite). In Mix I, complete stopping of balls in 
a shock occurs only in a minority of the shocks, but its 
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occurrence increases with polydispersity and x: In Mix 
III complete stopping of balls occurs in most shocks that 
have propagated for more than 20 — 40 cm and do not 
have another shock right in front of them and in Mix IV 
blocking occurs essentially in all shocks. 

To quantify this behavior further, we have measured 
(based on density space-time diagrams), for fixed values 
of (3 and D the fraction of time that the balls are stuck in 
such a stationary shock as a function of x (see Fig. |ll"|(b)). 
This data confirms what we already observed in Fig. llCj: 
the amount of blockage increases both with x and with 
the strength of the polydispersity. The increase with x 
can be understood simply from the observation that for 
strongly polydisperse mixtures shocks spread out dur- 
ing their lifespan, and since shocks travel upwards, the 
amount of blocked channel grows with x. 

One possible explanation we can find for the increase 
of this blocking with polydispersity is the occurrence of 
3D effects. One can imagine that in a shock wave, big- 
ger balls that are squeezed between small balls arc lifted 
from the support on which they roll. When such shock 
dissolves, small balls have to move over a finite distance 
before the bigger balls can start to roll, leading to a finite 
blocking time. 



Shock creation 



The polydispersity has two effects on the shock cre- 
ation. One could have anticipated that for stronger poly- 
dispersity the periodic packing sites become irrelevant, 
but as we see from Fig. [lj, what happens in addition is 
that all shocks are generated near the outflow of the fun- 
nel, an effect that can also be observed when comparing 
the density space-time diagram Fig. RG. 
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FIG. 12: Shock creation rates for monodisperse flow and poly- 
disperse mix I, II, III, IV. for f3 = 0.4°, D = 10 mm. 



2. Breakdown of scaling ofU(x) and v(x)D 

Using the RGE method we have studied to what ex- 
tend the scaling of U(x) and v(x)D with w(x)/D holds. 
In Fig. |l3|(a)-(b) we have plotted U(x) and v(x) for 
(3 = 0.4° and D = 10 mm for the four polydisperse mix- 
tures. For mix I, both U{x) and v(x) grow fairly lin- 
early with the usual packing site related variations ||. 
For mix II the behavior is similar but both U{x) and 
v(x) become flatter for large x. For the stronger poly- 
disperse mixtures III and IV strong deviations from the 
monodisperse or weakly polydisperse case are observed. 
The shock velocity U (x) grows rapidly for small x, then 
peaks and subsequently begins to drop for larger x. The 
shock frequencies v{x) in mix III and IV behave simi- 
lar to that of monodisperse/ weakly polydisperse flows at 
small x but become constant at larger x. This is fully 
consistent with the data in Fig. [l^ that shows that for 
strong polydisperse mixtures, all shocks are created near 
the outlet. 

Based on the similarity between weak polydisperse and 
monodisperse flows we plot U{x) vs. w(x)/D for mix I in 
Fig. |l3|(c). Apart from packing site variations the data 
collapse is reasonably good, so in this respect the weakly 
polydisperse flow behaves as a monodisperse flow. Linear 
fits of the data in Fig. ||(c) yield U(x) w 74[w(x)/D - 
0.18] cm/s. It is perhaps surprising that the mix I shock 
velocities are higher than those found for smooth wall 
monodisperse flows (see Section HA). 

The frequency deserves some more attention. When 
we plot v(x)D vs. w(x)/D in Fig. [b](d) the data col- 
lapse is not very convincing, especially for large x where 
the rescaled frequency is widely spread. This may be 
due to beginning effects of polydispersity (causing v to 
drop for high x as shown in Fig. [l^(b)) or partial failure 
of the RGE algorithm. This algorithm tends to discard 
weak shocks at high (3, x since the density contrast of 
the shocks becomes very small there (See §, @). bmce 
D = 7 mm data (highest shock frequency and lowest 




(b) P=0-4°, D=10mm 



B B { 

Mix I 
Mix II 
Mix III 

Mix IV 

_i i 



20 40 x(cm) 80 



20 40 x(cm) 80 



250 




2 w(x)/D 3 



2 w(x)/D 3 



FIG. 13: RGE based plots of U(x) (a,c,e) and v(x) (b,d,f): 

(a) shows U(x) vs. x at /3 = 0.4°, D = 10 mm for Mix I-IV. 

(b) shows u(x) vs. x at /3 = 0.4°, D = 10 mm for Mix I-IV. 

(c) shows U(x) vs. w(x)/D for Mix I. (d) shows v(x)D vs. 
w(x)/D for Mix I. (e) shows U(x) vs. w(x)/D for Mix III. 
(f) shows v{x)D vs. w(x)/D for Mix III. (In each of (c-f) 18 
data sets are displayed.) 



density contrast) show the biggest deviation the latter 
reason may be the most important. 

Despite the deviations for high x there seems to be 
a data collapse of v(x)D vs. w(x)/D which is su- 
perior to v(x) vs. w(x)/D: Together with t he re- 
analysis of the older data || discussed in Section II A 3 
this lead us to believe that the scaling of v(x)D with 
w(x)/D is a feature of both rough wall and smooth wall 
flows. Power Law fits of the data in Fig. |l^(d) yield 
v{x)D w l.6(l)(w(x)/D - l) - 7 ^) cm/s. That the shock 
frequencies of mix I are slightly lower than for monodis- 
perse flows is not surprising, since polydispersity gener- 
ally seems to make it harder for shock waves to form 
anywhere else than near the outflow of the funnel. 

For the U(x), v(x) statistics based on mix III density 
data shown in Fig. |l3](e)-(f) there is data collapse for 
neither U(x) nor i/(x)D vs. w{x)/D. All U(x) curves in 
Fig. [l3|(e) show the same pattern of rapid growth towards 
a maximum value after which a moderate decline sets in. 
We have found no clear pattern in the x and U values 
of the peaks, and there seems to be no w(x)/D scaling 



9 



involved. The v(x)D data for mix III shown in Fig. p"3|(f) 
displays growth at low x followed by a plateau, but there 
is no clear trend in this plateau value. It seems likely that 
both types of deviation from the monodisperse scaling 
relations are linked to the stagnant regions behind shocks 
and thus to the 3-dimensional packing effects discussed 
above. 

The U(x) and v(x) data for polydisperse mix II (not 
shown) are similar more to mix I, while the mix IV data 
(not shown) are similar to the mix III data. 

In conclusion, we find that for sufficiently strong poly- 
dispersity, the nature of the shock waves changes qualita- 
tively, and that scaling relations that hold for monodis- 
perse flows in either smooth or rough funnels break down. 



IV. DISCUSSION 

This work, in combination with earlier work on smooth 
wall/ monodisperse flows || |(|, leads to a number of 
conclusions about the effects of funnel geometry, wall 
roughness and inelastic dissipation. It is well-known that 
the dissipation occurring in ball-ball and ball-wall colli- 
sions is enhanced by the rolling nature of the ball motion 
|fL7| ; however, the precise value of the effective coefficient 
of restitution is presumably not of large importance for 
the phenomenology. In particular, the frequency and ve- 
locity scaling for monodisperse and weakly polydisperse 
systems are very similar (see Fig. 0), even though for 
weakly polydisperse flows dissipation seems enhanced. 
For strongly polydisperse flow, where dissipations seems 
very strong, we unfortunately loose the 2D nature of the 
experiment. 

In the case of smooth walls and monodisperse flows, 
packing effects become very important, and they tend to 
obscure scaling law. Possibly the simplest case is then 
the combination of rough walls and monodisperse flows 
(or, to a lesser degree, smooth walls and weakly poly- 
disperse flows). The rough walls have the additional ad- 
vantage that they make the system less sensitive to small 
perturbations (which are inevitably present for "smooth" 
walls) . 

The scaling laws for shock velocity and frequency, 
Eqs. are the main result of our work. Since these 
laws also approximately hold for smooth walls etc, their 
main origin must lie in the geometry of the experiment. 
That the shock velocity and frequency depend on x via 
w(x)/D is not surprising, since this is the most obvi- 
ous way in which x can be made non-dimensional (we 
do not expect the ball diameter to play an important 
role). The fact that v(x)D and not v scales is harder to 
understand, and may point at certain relevant velocity 
scales (note that both U and vD have the dimension of 
velocity). Two dimensional quantities that characterize 
the system are the effective ball-acceleration a e g (related 
to the inclination) and D, but these two quantities alone 
are not sufficient to provide for the correct scaling factors, 
since a velocity scale would be yJD a e ff, and a frequency 



scale would be ^/a^g/D, which both show scaling with 
D different from what is observed. 

We believe that the underlying reason for the occur- 
rence of the scaling laws is an important open question 
that deserves further study. Our data indicates that 
details of the ball-ball or ball-wall interactions are not 
important (although rough walls lead to more frequent 
shocks, they do not alter the nature of the scaling), which 
suggests that a relatively straight-forward model may 
capture the phenomenology here. Indeed, for pipeflow, 
where the geometrical cause of shock formation dissa- 
pears, leads to quite fragile behavior, which only in the 
case of rough walls seems to reproducable and similar to 
small funnel angle behavior. 

However, to make such a model one seems to need 
additional information about the relation between ball 
densities and velocities on the one hand, and shock fre- 
quencies and velocities on the other. 

Finally, shock waves have also been studied in some 
other geometries. Tsai et al. examined a two- 
dimensional rolling pipe flow of 3.2 mm steel spheres 
between rough walls. Partial blockage of the funnel out- 
let was used as an additional system parameter. These 
authors performed high resolution ball tracking mea- 
surements of small (12 cm) sections of the funnel and 
obtained detailed local measurements of v(x,t), p(x,t), 
T y (x,t), and v x (y). For an almost unblocked flow with 
D = 19 mm they found v ss 0.4 1/s and U « 14 cm/s, 
which is consistent with our data. 

Reydellet et al. J^] studied 1.5 mm metallic balls in a 
vertical pipe flow with rough walls. The balls were not 
rolling on a support, and the rolling of balls supposedly 
played less of a role than in our experiments. They qual- 
itatively observed the existence of upwards propagating 
shock waves, using a "double flash technique" to study 
local ball velocities and made measurements similar to 
those in our Fig. |9[ Note that their funnel was not con- 
tinuously refilled from a reservoir and thus their studies 
were limited to the transient behavior. All our experi- 
ments are preceded by a > 30 sec. flow in order to avoid 
such transient behaviors. 

Finally, Le Pennec et al. Q studied a 2-dimensional 
rolling flow of 1 mm glass balls in flow geometries with 
very large /3 (mostly 30°), large D (typically 10-120 ball 
diameters) at various flow plane inclinations. Due to the 
different geometry it is hard to make any direct compari- 
son, although these authors did measure shock velocities 
and frequencies. 



Summary 

The packing effects that we observed in earlier work 
on monodisperse, smooth walled flow ||] can be sup- 
pressed by either making the sidewalls rough or using 
polydisperse mixtures of balls. For rough walls, we find 
that there are more shock waves than for smooth walls as 
a result of the greater disturbances in this flow. In earlier 
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work H on smooth wall flows we found that the average 
shock speed U(x) scales with w{x)/D. We have found 
a similar scaling for rough walls and observed that also 
v(x)D scales with w(x)/D for both rough and smooth 
walls 
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By using ball tracking methods we shed light 
on the shear flow properties for rough walled flows. 

In polydisperse flows between smooth walls we find 
that our scaling relations persist for very weak polydis- 
persity, but breaks down for stronger polydispersity. This 
breakdown is most likely caused by weak 3-dimensional 
effects in the packing of shock waves which apparently 
lead to partial blocking of the funnel. By using ball track- 
ing methods we shed light on the shock structure and on 
the extended stationary shock packings found in strong 
polydisperse flows. 
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APPENDIX: COMPARISON WITH TRAFFIC 
FLOW 

To compare our flows to traffic flows, small angle, or 
preferably pipe flows should be considered. In earlier 
work we made comparisons between smooth wall funnel 
flows with nonzero (3 and traffic flows |(| . In smooth wall 
pipe flows, shock waves are extremely fragile, and such 
a comparison is not meaningful. For the rough wall pipe 
flow discussed here, however, such a comparison can be 
made in principle. In addition, we will discuss our results 
for rough wall funnel flows in terms of traffic flows. 

It is commonly assumed that there are three flow types 
in traffic flow, namely, uncongested flow (a steady flow of 
vehicles at low to moderate densities, queue flow (a slow 
flow of near maximum density), and queue discharge (a 
flow of vehicles accelerating out of a queue flow) M . 
Traffic data are usually represented as v(q) (speed/flow 
relation) or as q(p) (fundamental diagram), where q = pv 
is the flow rate. 

In Fig. |l4|(a)-(b) we show grey scale histograms of 
(q, v) and (p, q) for a rough wall pipe flow (j3 = 0°, 
D = 15 mm). Note that we use the relative density p 
instead of p. The average values, corresponding to v(q) 
and q(p), are shown as solid lines. In Fig. |lj(a), we 
observe regions corresponding to queue flow and queue 
discharge, but there is no region corresponding to uncon- 
gested flow. As shown in Fig. |lj(b), the flow rate q(p) 
has a parabolic shape with a peak value around p = 0.6. 

For larger values of D (not shown), the flow rates ex- 
hibit similar parabolic behavior, albeit with larger typical 
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FIG. 14: Fundamental diagrams in rough wall flows, (a-b) 
Pipe flow studied in the interval 37 < x < 74 cm for = 
0.0° and D — 15 mm. (c-d) Intermittent flow studied in the 
interval < x < 37 cm for (3 = 0.4° and D — 20 mm. 



values ofw and q, and with the q(p) "parabola" skewed to- 
wards lower p (p = 0.4-0.5 for D = 20 mm, p = 0.2-0.4 
for D = 25 mm). 

Figs. |l4|c-d show the corresponding diagrams for small 
angle, rough wall funnel flow at (3 = 0.4°, D = 20 mm. 
We find queue flow and queue discharge but no uncon- 
gested flow, similar to what we found in earlier work on 
smooth wall funnel flow Q. In Fig. |l4| we observe a 
parabolic shaped q(p). Similar plots for higher f3 (not 
shown) exhibit "parabolas" skewed towards higher val- 
ues of p. 



APPENDIX: VELOCITY AND FREQUENCY 
DETERMINATION FOR PIPEFLOW 



The method to obtain shock velocities for pipeflow 
data is illustrated in Fig. [l5| First of all, to highlight 
the shocks that occur in the density plots, some smooth- 
ing and differentiation is applied yielding a field s(x, t) 
as shown in Fig. [l5](a). Clearly, shocks are now visi- 
ble as bright streaks (high values of s) in a fairly even 
background. The temporal averages of the spatiotempo- 
ral correlation function C(x 1 Ax 7 At) := J dt s(cc, t)s(x + 
Ax, t+At), two examples of which are shown in Fig.[lq(b- 
c) for x=130 and 170 cm, clearly show a dominant di- 
rection in space-time that can be associated with the 
local dominant velocity of the shocks: this method of 
shock velocity determination works over the whole range 
of parameters considered. In comparison to the RGE 
method, this method is more local, and often gives a 
slightly smaller estimate for the velocities (order of 10-20 
%). This is presumably due to the fact that shocks have 
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a tendency to perform intermittent jumps forward, mak- 
ing their dominant local velocity smaller than the typical 
velocity obtained over longer timescales. Some effect of 
these jumps shows up in the correlation functions, where 
the bright streak tends to bend for larger correlation dis- 
tances and time intervals. 
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FIG. 15: Illustration of methods used to obtain shock veloci- 
ties and frequencies for pipe flow for the case of D — 20 mm. 
(a) Modified density plot, (b-c) space-time correlation func- 
tions based on data shown in (a), for (b) x = 130 cm and 
(c) x = 170 cm. Note that the scales of (a) and (b,c) are 
different, so that the dominant angles appear different, (d) 
Thresholded data used for determination of shock frequency. 



To obtain the frequency from s(x,t) is fairly straight- 
forward: after choosing an appropriate threshold, one 
obtains two-color images like shown in Fig. |ll](d). A 
simple algorithm suffices to count, for a fixed value of x, 
the number of shock waves that occur. 
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The shock detection algorithm RGE (.Relative density 
contrast Gradient Edge detection) used on p-data is a re- 
finement of the GE- method described in || . It is used to 
measure U(x) and v(x). We use two versions: RGE4 with 
an a>res. of 25 cm and RGE8 with an a;-res. of 12.5 cm. 
The RGE algorithm gives more reliable measurements 
of average shock frequency v{x) but may still underes- 
timate it under conditions with high v and very low p 
contrast. The RGE method does not detect shocks with 
U < 15 — 20 cm/s as reliably. This is only a real problem 
in rough walled pipe flows. 

This was not shown in Q], but with improved data anal- 
ysis jl4 of the old data we find a decent data collapse 
when p otting v(x)D vs. w(x)/D. 

The Shock Creation Rates shown here have been hand 
counted, and every shock that propagates upwards at 
least 5 cm (no matter its strength or continued existence) 
is counted. The shock wave creation position statistics 
shown in Fig. 18 in |f| were based on a image processing 
method (DT - no longer used), that only counted shocks 
of a certain strength. Consequently this method showed a 
decreasing shock creation rate with increasing w(x) (for 
smooth walls). 
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